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Abstract 

Dromions are exponentially localised coherent structures supported by nonlinear integrable 
evolution equations in two spatial dimensions. In the study of initial-value problems on the plane, 
such solutions occur only if one imposes nontrivial boundary conditions at infinity, a situation of 
dubious physical significance. However it is established here that dromions appear naturally in the 
study of boundary-value problems. In particular, it is shown that the long time asymptotics of the 
solution of the Davey-Stewartson I equation in the quarter plane with arbitrary initial conditions 
and with zero Dirichlet boundary conditions is dominated by dromions. The case of non-zero 
Dirichlet boundary conditions is also discussed. 
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1 Introduction 



We consider the Davey-Stewartson (DS) I equation [1] 

m + ^(fe + Qyy) - {<Px + \qf)q = 0, 

<Pxx - fyy + 2|gpX = 0. (1.1) 

In the context of water waves this equation is the shallow water limit of the Benney-Roskes [2] 
equation in the case of dominant surface tension. In this case q{x, y, t) is the amplitude of a surface 
wave packet and ip is the velocity potential of the associated mean flow. In general the DSI equation 
provides a two-dimensional generalization of the nonlinear Schrodinger equation and can be derived 
from general asymptotic considerations [3]. 

We introduce characteristic coordinates and we also replace the second order equation (1.1b) 
by two first order equations: Let ^, ry, Ui{^, rj, t), U2{i, i], t), be defined by 

1 1 
^ = x + y, r] = x-y, Ui = -tpri - ^\q\'\ U2 = -tp^ - -\q\'^ . (1.2) 

Then the DSI equation becomes 

iqt + q^^ + qrm + (^i + U2)q = 0, 

Ui, = \\q\l U2^ = \\q\l (1.3) 

Indeed, writing equation (1.1 a) in characteristic coordinates and using the definitions of Ui,U2, 
we find equation (1.3a). Also the definition of C/i implies C^i^ = —^t]^ — klf/2; writing equation 
(1.1b) in characteristic coordinates and using this equation to replace y?^^, we find equation (1.3b). 
Similarly for equation (1.3c). 

Formulation of the Problem. Let the complex-valued function g(^, ?7,t) and the real-valued 
functions Ui{^,r),t), U2{^,rf,t), satisfy equations (1.3) in the domain 

0<^<oo, 0<r?<oo, t>0, (1.4) 
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with the following initial and boundary conditions: 



q{C,V,0) = Qo{^,v), 0<(<oo, 0<r?<oo, 



gr(0, r?, t) = 0, Ui{0, r?, t) = «i(r?, t), < ry < oo, i > 0, 



g(e,0,t) = 0, C/2(C,0,t) =«2(e,i)> 0<C<oo, t>0, 



(1.5) 



where the functions qo{C,r]), 1x1(77, t), U2(C,i) have sufficient smoothness and they also decay as 
^ ^ 00 and 77 ^ 00. 

Notations. 

• Bar denotes complex conjugation. 

• Mil, M12, M21, M22 denote the (11), (12), (21), (22) entries of the 2 x 2 matrix M. 

• Md^Mq denote the diagonal and the off-diagonal parts of the 2x2 matrix M. 

Theorem 1.1 Given qo{£,,ri) define the vector {Mi(^,r],k), M2{^,r],k))^ by 

Mij - ikMi = ~qoM2, = ^QqMi, 0<^<oo, 0<r/<oo, Imk < 0, 



^— >oo 



limMi = 0, M2(C,0,A;) = 1. 



(1.6) 



Given qo and M2 define SQ{k,l) by 




ImA; < 0. 



(1.7) 



Given Sq, ui, U2, define S{^,r),t) by 



iSt + S^^ + Snr^ + {ui{r],t) + U2{^,t))S = 0, < ^ < 00, 0<r/<oo, t > 0, 




< ^ < 00, <rj < 00, 



S{^,0,t) = 0, 0<^<cx), t>0, 
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S{0, r),t) = 0, < r/ < oo, t > 0. 



(1.8) 



Given S define the 2x2 matrices M~^{^,r],t,k) and M {^,r],t,k) as the solution of the following 
non-local Riemann-Hilbert problem: 

• M+ and M~ are analytic for ImA; > and ImA; < respectively. 

• M± = / + 0(i), A;^oo, Im/cT^O. 



For real k, M+ and M satisfy the following jump condition, 



I 



i^,r],t,k) - 



/oo 
dlS{l,k,t)t 
-oo 



-il^—ikr] 



12 



V ^22 



/ 



i^,rj,t,k) = - dlS{Kl,t)e''^+'^^ 



M- 



11 



V ^2l 



where 



1 ;>oo roo ^ 

S{k,l,t) = ^J^ d^dvSi^,ri,t)e-'''^-''\ 



i^,r],t,l), (1.9) 



(1.10) 



This Riemann-Hilbert problem has a unique global solution. 
Define q{^,v,t) by 



q = 2i lim {kMi2). 

k—KX> 



(1.11) 



Then q satisfies (1.3), (1.5). 



Remark 1.1 Although the evolution of the scattering data S{k,l,t) is in principle determined 
by equations (1.8), the relevant time dependence is complicated. In turn, this makes it difficult 
to determine the long time behavior of the solution M of the RH problem (1.9). This difficulty 
can be bypassed by formulating an inverse problem for S{k,l,t). Since equation (1.8a) is precisely 
the one studied in [4] , the relevant analysis is identical with the one presented in [4] : (a) If ui and 
U2 are time-independent then the analysis of (1.8) is intimately related with the analysis of the 
time-independent Schrodinger equation 
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where u is either ui or U2- After a long time the solution of (1.8) is dominated by the discrete 
spectrum of ui and of U2- If ui and U2 have A^i and N2 discrete eigenvalues respectively then the long 
time asymptotics of q is given by an (iVi, N2)-breather solution, (b) If ui and U2 are time-dependent 
then the analysis of (1.8) is related with the analysis of the time-dependent Schrodinger equation. 
Assuming a certain completeness relation for the eigenfunctions of the time-dependent Schrodinger 
equation, the long time asymptotics of the solution q is dominated again by the associated discrete 
spectrum; in this case the solution is dominated by an {Ni, N2)-dromion solution. 

Remark 1.2 The explicit form of the (iVi, A'^2)-dromion solution can be found in [4], see also 
[5], [6]. Dromion solutions are exponentially decaying in both ^ and 77. However, in contrast to the 
one-dimensional solitons, these solutions do not preserve their form upon interaction and therefore 
can exchange energy. These coherent strucures can be driven everywhere in the quarter-plane 
C ^ 0)^ ^ 0) by choosing a suitable motion for ui{ri,t),U2{£,,t). 

Remark 1.3 A characteristic feature of boundary- value problems for integrable evolution PDEs 
in one spatial dimension is that the (spectral) functions defining the associated Riemann-Hilbert 
problems cannot in general be expressed explicitly in terms of the given boundary conditions [7]. 
Explicit formulae can be obtained only for a particular class of boundary conditions which are 
referred to in [7] as linearisable. The situation is similar for boundary-value problems in two spatial 
dimensions. We emphasize that equations (1.8) are explicitly defined in terms of the given initial 
conditions, thus the homogeneous Dirichlet problem (1.3)-(1.5) belongs to the class of linearisable 
conditions. It is shown below that for non-homogeneous Dirichlet boundary conditions, the equa- 
tions defining S, in addition to the given Dirichlet boundary conditions, also involve the unknown 
Neumann boundary values. 

Proposition 1.1 Let the complex-valued function q{i,ri,t) satisfy a non-homogeneous Dirichlet 
boundary-value problem in the quarter plane, namely equations (1.3)-(1.5), where q(0,r],t) = 
and ^(^,0,^) =0 are replaced by 

q{0,ri,t)=go{ri,t), g(e, 0, t) = /o(C, t), (1.12) 



5 



and the functions go, /o have sufficient smoothness and they also decay as 77 — > oo, ^ — > oo. Assume 
that there exists a global solution. 

This solution can be expressed by equation (1.11) through the solution of the Riemann-Hilbert 
problem defined in Theorem 1.1. This problem is uniquely defined in terms of the function S(^, r], t) 
which satisfy the following: 

• >S' solves the linear evolution PDE 

iSt + S^^ + Sr,r, + {ui + U2)S+ c?|Fi (e, |, 7?, + P dfiF2{r),fi,t)Si^,fi,t) = 0, (1.13) 

Jo Jo 

where 

m,lt) = i[h{i,t)fo{c,t) - fo{it)h{c,t)] - y{^,t)f{it)jide\f{e,t)\'', 

F2{ri,fi,t) = l[gi{fi,t)go{v,t) - go{fi,t)gi{ri,t)] - ^giri,t)g{fi,t) dri'\g{r)' ,t)\'', 

and gi{rj, t), fi{^, t) denote the Neumann boundary values q^{0, 77, t), qr]{^, 0, t); 

• S satisfies the initial and boundary conditions 

S{^,V,0) = Soi^,v), S(C,0,t)=7r/o(C,t), S{0,ri,t) = 7rgoiri,t), (1.14) 
where So{^,i]) is defined in terms of qo{^,v) by the rhs of equation (1.7). 

Remark 1.4 In the case of boundary-value problems for evolution PDEs in one spatial dimension, 
the unknown boundary values can be characterised in terms of the given boundary conditions 
through the analysis of certain global relations [7], [8]. Such global relations exist in two spatial 
dimensions and can also be used for the characterisation of the unknown boundary values. This 
analysis, which is rather complicated, will be presented elsewhere. 

Organisation of the paper. The DSI equation admits a Lax pair formulation. The t-independent 

part of the Lax pair (sec equations (4.2)) is analysed in section 2. The specific form of the t-part 
of the Lax pair (see equations (3.4)) depends on the specific form of the boundary conditions. The 
t-paxt of the Lax pair for the homogeneous and non-homogeneous Dirichlet cases is discussed in 
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sections 3 and 5 respectively. The proof of Theorem 1.1 is presented in section 4. Section 6 contains 
further discussion. 



2 The independent part of the Lax Pair 

Throughout this section we suppress the t-dependence. The t-independent part of the Lax pair is 
given by equation (4.2). Analysing this equation in characteristic coordinates we find the following: 

Proposition 2.1 Let the vectors 



M+ M+ Mfi 



(2.1) 



(2.2) 



which are functions of ^, 77, k, be defined by 

M^i = -\ /r drj'e^''(^'-^kM^„ M+ = 1 + i /J' dv'qMt2, 
Mn = 1-5/0 d^qM^,, Mf^ = ^ de e'''(^-^") qM^^, 

M21 = I IS dq'qe^'^^^'-^^Mri, M^^ = 1 + l /J' dv'qM^^, 

where q{^,rf) has sufficient smoothness and decay. Then: 

1. The first two vectors in (2.1) are analytic in k for ImA; > 0, while the last two vectors in 
(2.1) arc analytic for ImA; < 0. 

2. For real k, the vectors (2.1) satisfy the relationships 
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/ 








V 


j 




V 


M22 1 



dlSik,l)e^'^+^'^n \{^,V,l), (2.3) 



21 



where S{k,l) is defined by 

1 roo roo 

Sik, = — / d^dm{i. V)M22 (e, T?, fc)e-^'=«-^''', ImA; < 0, ImZ < 0. (2.4) 
47r Jo Jo 

3. The vectors (2.1) have the following behaviour for large k, 

Mti = i + o(^^y M2± = l + oQ), M±2 = oQ), Mti = o(^^y k^oo, Imfe ^ 0. 

(2.5) 

Proof 

1. The vector {M^^^, M2i)'^ satisfies a system of linear Volterra integral equations with kernel 
analytic in k for Irak > (since rj' > ij). Thus this vector is analytic for Irak > 0. Similarly for the 

other vectors. 

2. Let the matrices ^~^{^,r],k) and "^~{^,r],k) be defined by 



'■22'' 



(2.6)= 



Then 



■22' 



*2"l = 5 /O ^^V^-^Fl, *2"2 = e-^'^^ + \ IS dv'm2- 



(2.7) 
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Subtracting the equations defining the vectors (^i2)*22)'^ ^^'^ (^125^22)"^ ^^'^ 



1 POO 1 
*^2 - *r2 = "2 X ^^'^*22 - 2 X ^^'^^^^2 - *2"2), 



^t2 - ^2"2 = 2 /J ^'?'^(*^2 - *r2)- (2. 



Using the definition of S{k, I) it follows that the first term of the rhs of equation (2.8a) equals 



/oo 
dlS{k,l)( 
-00 



Mr) 



Comparing equations (2.8) with the equation satisfied by {^11, *2i)^ equation (2.3b) writ- 

ten in terms of the functions instead of the functions M^. 

Similarly subtracting the equations defining the vectors (*ii,*2i)'^ and (^n'^21)'^ 



equation (2.3a), where instead of —S{l,k) we have the function 



"I POO POO 

Tik,l) = -—J^ io didm{i,r^)M+{i,r^,ky 



(2.9) 



We will now show that 



T{k,l) = -S{l,k). 



(2.10) 



Indeed, the vectors (^ii,*2i)'^ (^125^22)'^ satisfy the equations 



12. 



22 



*2l5 = 2^*11' 



Hence 



22: 



n 2 



o9*12- 



(*+(/c)*r2(0) =-(*2+i(^)*2"2 



where for convenience of notation we have suppressed the (^, rj) dependence. Integrating this 
equation we find 

Jo L 



^=00 
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ri=0 



■n=o 



Using equations (2.7) to compute the boundary values appearing above, for example 



12 



^=00 
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= 0, ^r+(fe) e^fc", *22(0 

•q=oo q=U 



ri=Q 



we find 



Jo Jo Jo Jo 

The Ihs of this equation is the complex conjugate of S{l,k), while the rhs equals —T{k,l). 
3. Equations (2.2) and integration by parts, imply equations (2.5). 



QED 



3 The t-part of the Lax Pair for the homogenous Dirichlet case 

In what follows we first derive the t-part of the Lax pair for the vector {"^12^ ^22)^- 

Proposition 3.1 Assume that there exists a function q{C,r),t) with sufficient smoothness and 
decay which satisfies equations (1.3)-(1.5). Let the vector * = (^'i,^'2)'^ satisfy 

1 r°° 

^i{^,ri,t,k) = -y^ d^'qi^',ri,t)^2{^',ri,t,k), 

^2i(,V,t,k)=e-"'^ + l r dri'q{^,v',t)^i{^,v',t,k), Imfc < 0. (3.1) 

2 Jo 

Then the function tpi(ri, t, k) defined by 

Mv,t,k) = ^i{0,r],t,k), 

solves 



1 f'°° 

^V'lt+V'i,, +«i(??,i)V'i-(^'V'i + o9«(0'^'*))-^ / dlj{k-l,t)Mri,t,l) = 0, (3.2) 

J —00 



10 



where 



7 /'OO 

^ik,t) = -J^ d^U2i^,t)e-"^^. (3.3) 



Proof We will first show that if ^2)^ satisfies equations (3.1) and q{^, r], t) satisfies equations 
(1.3)-(1.5),then, 



= i{d^ - 5^)2^'i + iq{d^ - 5^)5'2 + iC/i^'i - + vi, 



(3.4) 



*2t = -i{d^ - dr^f^2 - iq{d^ - ^j?)*! - iC^2*2 - + V2, 
where the vector v = {vi,V2)'^ is given by 

/OO 
dh{k-l,t)^{^,ri,t,l). (3.5) 
-00 

Indeed, the vector * = (*i, ^'2)^ satisfies 

= -^5*2, *2^ = (3.6) 

Suppose that ^ satisfies equations similar with (3.4) but with v = 0. It is straightforward to 
verify that the compatibility condition of equations (3.4) with v = 0, and of equations (3.6) yields 
equations (1.3). Actually, if equations (1.3) are valid then equations (3.4) and (3.6) are compatible 
for any vector v given by 

/OO 
dlT{k,l,t)^{^,r],t,l). 
-00 

The precise form of F depends on the boundary conditions that ^' and q satisfy. Indeed, ^'i can 
be writing in the form 

/■OO roo 

*i = -y or = -y dei^i,)^'. 
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By replacing in the above equation hy A + vi, where A is defined in the rhs of equation (3.4a), 
we find 



vi = — Um 
5— >oo 



Using 



hm q = hm ^'i = 0, 

^— >oo ^— »oo 



the term in the bracket vanishes. Also 



/oo roo / 

cZ/r*i^ = / dir[- 
-oo J —oo \ 



Hence 



Similarly starting with 



1 
2 



(3.7a) 



Jo ^ Jo 



we find 



Using 



V2 



1 If 

i(5£-9„)2^'2+i^(9£-9„)*i + if/2^'2 + %*i +7:/ drj'qvi. 

J »7=0 2 Jo 



it follows that 



1 f^i 

V2 = -ik'^e-'''^ + iu2{i, t)e-'^^ + - / dr]'qvi. 

2. Jo 

Equations (3.7) imply that the term — i/c^ exp[— i/c^] gives rise to —ik^'^, while the term 



(3.76) 



/oo 
dl^{k-l,ty 
-00 



-Hi 



(3.8) 



gives rise to the term J_^dlj{k — l,t)'^{l). Also equation (3.8) implies that 7(A;,t) is given by 
equation (3.3). 
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Equation (3.2) follows from the evaluation of equation (3.4a) at ^ = 0. In this respect we note: 

1 



i=o 2 



since q{0,ri,t) = q^{^,rj,t) = 0. Also 



^1, 



-\ [q^^2 + g«'2d£=o = -^9^(0, ??, t), 



since *2(0, 77, t, /c) = 1. 



QED 



Using equation (3.2), it is now straightforward to obtain the time evolution of S{k,l,t): 

Proposition 3.2 Let S{k, I, t) be defined by equation (2.4), where M2'^(^, r], t, k) = ^'2 exp [— ifc^], 
^'2 is defined in terms of q{$,,r],t) in equations (3.1), and q evolves in time according to equations 
(1.3)-(1.5). Let 

/oo poo 
/ dkdle'''^+'^'^S{k,l,t). (3.9) 
-oo J —oo 

Then S satisfies equations (1.8). 

Proof. Using (M{^,M2^)^ = e-^'=^(*i, ^2)^, it follows that the vector = (*i,*2)^ satisfies 
equations (3.1). Also 

1 roo poo 1 l-OO 

S{k,l,t) = —j^ io ^^'^W(?'^)^"'"'*2(?,?7,fc) = ^y^ dr?e-''^*i(0,r?,t,fc), (3.10) 

where the first equality uses the definition of S (see equation (2.4)) and the second equality uses 
equation (3.1a). Hence using ^'i(0, 77, t, k) = ipi{r], t, k) and replacing in equation (3.9) S by the rhs 
of equation (3.10) we find 

S{^,ri,t)= dke"'^Mri,t,k). (3.11) 
J —00 
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In order to find the time evolution of S we need to take the A;-Fourier transform of equation (3.2). 
In this respect we note that if q{0,'rj,t) = 0, equations (3.1) imply 



A; ^ oo. 



(3.12) 



Using this estimate, equation (3.2) implies equation (1.8a). In this respect we note that the k- 
Fourier transform of k^tpi + ^^'^(0, 77, t) equals S^^. Indeed, if s{x) denotes the Fourier transform of 



/oo 
e'^''ij{k)dk, X > 0, 
-00 



and if il){k) for large A; is given by 



then 



^= e'^''[a-k^'il:{k)\dk. 
In order to derive this result we note that the definition of s(x) can be rewritten as 



/ex 
-c 



Akx 



i^{k) + 



a 



1 + A;2 



/oo 



dk; 



the second integral above can be computed explicitly and equals ire ^ . Thus 



dx^ 



-ikx 



ak'^ 



l + P 



dk + aire ^; 



simplifying the rhs of this equation we find (3.15). 

Since V'i(0,i,/c) = 5'i(0,0,i,fc) = 0, equation (3.11) implies 5(^,0,*) = 0. Also 



/ex. 
-c 



dkilji{r], t, k) = 0, 



since ■^i is analytic for ImA; < and is of 0{X) as k ^ 00 (see equation (3.12)). 



(3.13) 



(3.14) 



(3.15) 



QED 
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4 Proof of Theorem 1.1 



In both sections 2 and 3 it was assumed that rj, t) exists, it is smooth and has sufficient decay. 
Then, under this assumption, it was shown in section 2 that q can be obtained through the solution 
of a RH problem uniquely defined in terms of S{k,l,t). Furthermore, it was shown in section 3 
that S satisfies equations (1.8). We will now prove that q{^,'r],t) can be constructed through the 
solution of the RH problem (1.9) without the a priori assumption of existence. Furthermore, this 
result is valid without the need for the norms of qo{^,t),ui{rj,t),U2{^,t) to be small. 

The function So{k,l) is defined through the solution (Mi(^, 77), M2(^, r/))-^ of the linear integral 
equations (1.6). These equations are of Volterra type, and thus if qo e Li, are always solvable. 

The Fourier transform S{^, r], t) of the function S{k, I, t) satisfies the linear PDE (1.8a) with the 
initial condition S'o(C, rj) and with homogeneous Dirichlet boundary conditions. Thus if ui{'q, t),U2{^, t) 
have sufficient smoothness and decay, S{k,l,t) is well defined. 

If the functions S, Sk, Si G L2 for fixed t, the RH problem (1.9) has a global solution provided 
that its homogeneous version has only the trivial solution. This is indeed the case [4]: Let satisfy 
equations (1.9) but with /j, = 0{l/k) as k ^ 00. The jump conditions imply 



where for convenience of notation we have suppressed the ^, rj, t dependence of jU, S. Multiplying 
equation (4.1a) by ^^^(k), equation (4.1b) by ^22(^)^ integrating over k, and adding the resulting 
equations we find an equation whose Ihs equals 



while the rhs is imaginary (in this derivation we have used that the integrals of A*ii(fc)A'ii(fc) and 
of iJ.22(k) Hi2i^) vanish since these products are analytic for ImA; > and ImA; < respectively, and 



For the proof of the fact that q{^,r],0) = qo{^,v) we note that q{^,rj,0) is defined through the 
RH problem (1.9) but with So{k, I) instead of S{k, 1, 0). By repeating the analysis of section 2 where 



/^ii(^) -/^ii(^) = 

/^12(^) -/^r2(^) = 



JZ,dlS{l,k)e-^^i-^'^fit2il), 
!^^dlS{k,l)e^'^+'^^H^^{l), 



(4.1) 




are of 0{l/k^) as k ^ 00). Hence, iJ.^ = fii2 = 0- 
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q{^,r),t) is now replaced by the given function q'o(C)'?)) it follows that qo{^,rj) is characterised by 
precisely the same RH problem as the problem characterising q{^, 77, 0). (We note that (Mi, M2)'^ = 
{M12, M22)'^)- Hence the unique solution of this RH problem implies g(C,??, 0) = qoiC^i])- 

We now show that if q{^,r),t) is defined through the solution of the RH problem (1.9) then q 
satisfies the DSI equation. This proof is based on the extended version of the dressing method 
presented in [14]. The dressing method can be used to show that if M solves an appropriate RH 
problem then M also solves both parts of the associated Lax pair. Hence, using the compatibility 
condition of this pair, q solves the relevant nonlinear PDE. 

We first show that M solves the t-independent part of the Lax pair 

M^ + a3My-ik[a3,M] + QM = 0, (4.2) 

where 

( » A 

i = x + y, rj = x-y, (73 = diag(l, -1), Q= , (4.3) 

\-q ; 

and [, ] denotes the usual matrix commutator. Indeed, the jump condition of the RH problem can 
be written in the form 

/oo 
dlM{x,y,t,l)F{x,y,t,k,l), (4.4) 
-00 

where M denotes the matrix with first, second column vectors (Mj^,M^)'^, (M{2,M^)^ respec- 
tively, and F denotes the off-diagonal matrix with 

Fi2{x,y,t,k,l) = -S{k,l,t)e^^^''+y^+^^^^-y\ 

(4-5) 

F2i{x,y,t,k,l) = Fi2{x,y,t,l,k). 

Using 

Fi2x = i{k + l)Fi2, F2i^ = -i{k + l)F2i, 
Fi2y = i{k-l)Fu, F2iy = i{k - l)F2i, 
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and writing these equations in a matrix form we find 



= ilasF - ikFas, Fy = i{k - l)F. (4.6) 
Let the operators Dx,Dy be defined by 

D^M = Mx + ikMa3, DyM = My -ikM. (4.7) 
These definitions and equations (4.6) imply that 

dlM{l)F{k, l)= dl (d^M{1)) F{k, I), 

-oo J—oo 

/oo /* oo 

dlM{l)F{k, l)= dl (DyM{l)) F{k, I), 
-oo J—oo ^ ' 

where for simphcity of notation we have suppressed the x, y, t dependence. Hence, both D^M and 
DyM satisfy the same jump condition as M. Using the fact that QM also satisfies the same jump 
condition, it follows that the expression 

D^M + asDyM + QM (4.8) 

satisfies the same jump condition as M. Using 

^ = ^ + -k^ fe^ + ^lpj' (^-^^ 

it follows that the large k behaviour of the expression (4.8) is given by 

(g-i[a3,MW])+oQ). 

Hence if we define Q by 

Q{x,y,t)=i[u3,M^% (4.10) 

the expression (4.8) satisfies the homogeneous version of the RH problem, therefore it vanishes and 
M satisfies equation (4.2). We note that Q is off diagonal; furthermore equation (4.5b) implies 
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certain symmetry relations for M, which in turn imphes that {Q)2i = —{Q)i2- 1^ summary, if M 
satisfies the RH problem (1.9) and if q is defined through M by equation (1.11), then M satisfies 
the t-independent part of the Lax pair. 

We now show that M also satisfies the t-pait of the Lax pair. For this purpose, rather than 
assuming that S satisfies equation (1.8a) we first assume that S{k,l,t) satisfies the equation 

/oo 
du{-i)ji{u - I, t)S{k, u, t) + Z72(i/, t)S{k - u, I, t), (4.11) 
-oo 

where the functions 71 , 72 are defined by 

j fco j roc 

M^,t) = — / d^?nl(77,^)e-^ j2(.iy,t) = -- / den2(^,i)e--«. (4.12) 

^TT J-00 ^TT J-00 

These equations imply that F satisfies the evolution equation 

Ft=ik'^Fa^-ifa^F+ H du(r{x,y,t,u-l)F{x,y,t,k,u)-F{x,y,t,k-u,l)r{x,y,t,u)) , (4.13) 

J —00 

where T is defined by 

r(x,y,t,z.) = e-^diag(e-*^^7i(^^),e'''"72(z^)) • (4.14) 

Indeed, using 

Fi2(x,y,i,fc,0 = -5(fc,^,^)e*'=«+*'^ Tnix,y,t,iy) = e-''''^ji{u,t), T22{x,y,t,u)=e"'^j2{i^,t), 

the (12) component of equation (4.13) yields equation (4.11). 

The time dependence of F suggests the introduction of the operator Dt defined by 

/oo 
dvM{x,y,t,k-u)T{x,y,t,v). (4.15) 
-00 

Indeed it can be shown that 

/oo /'OO 
dlM{l)F{k,l) = dl{DtM{l))F{k,l). 
-00 J— 00 
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This equation is a direct consequence of the definition of Dt and of equation (4.13) provided that 



/OO POO POO POO 

dl dvM{l)T{v-l)F{k,u)= dl dvM{l - v)T{v)F{k,l)- 
-OO J — OO J — OO J — OO 

this equality is established by replacing in the rhs I — v hy v and then exchanging I and v in the 
resulting equation. 

Since DtM satisfies the same jump condition as M, following the same logic as in the derivation 
of equation (4.2) it follows that 

-DtM + ia^Dy^M + BDyM + AM = {), (4.16) 

where the matrices A{x^ y, t) and B(x, y, t) will be chosen by the requirement that the Ihs of equation 
(4.16) is of 0{l/k) as k ^ OO. In this respect we first introduce the following notations for the 
0{l/k) term of M, see equation (4.9), 



^ a b 



c d 



(4.17) 



M^^\x,y,t) = 
Equation (4.10) yields 

The 0{l/k) term of equation (4.2) implies 

MW + agM^^) - i[a3, M^^)] + QM^^ = 0. (4.19) 
The diagonal and off-diagonal parts of this equation yield 

9^Mg^ + (T^dyM^j^^ + QM^') = 0, (4.20) 

d^M'^O^ + a^dyM'^^^ - i[(73, Mg'^] + QMg^ = 0. (4.21) 



Equation (4.20) implies 
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Using the definitions of Dt and Dy, equation (4.16) becomes 

Mt = -ik'^[a3, M] - ikBM + 2ka2,My + AM + BMy + ia^Myy 

/oo 
duM {x, y,t,k- iy)T{x, y, t, u). (4.23) 
-oo 

The 0{k) term of this equation yields 

S = -[a3,M«], 

thus (comparing with equation (4.10)) 

B = iQ, (4.24) 

The 0(1) term of equation (4.23) yields 

/oo 
dvV{x,y,t,v)=0. (4.25) 
-oo 

Using equation (4.19) to replace the first two terms of this equation we find 

/oo 
dvT{x,y,t,v) = Q. 
-oo 

Solving this equation for A and using the definition of T (equations (4.12)) we find 

All = 2ar^ + iui, Ai2 = 2b^, A21 = 2c^, A22 = 2d^ - iu2. (4.26) 

In summary, if the time evolution of S{k,l,t) is given by equation (4.11), where 71 and 72 are 
defined by equations (4.12), then the time evolution of M is given by (4.23) where B = iQ and the 
components of the matrix A are given by equations (4.26). 
Letting 

M = ^diag (e-''''^, e'^^) , (4.27) 
and using that the y-derivative of the above diagonal matrix equals ik times this matrix, it follows 
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that the time evolution of * is given by 

-^t = i(^^dl^ + iQdy^ + + ik'^^(j^- j dz/*(a;,y,t, A; - i/)diag(7i(i/,t),72(i^,i)). (4.28) 
With a particular choice of the constants of integration, equations (4.22) yield 

Then the first and the fourth of equations (4.26) become 

All = ^{lj^ + MV, t)j , A22 = -i Q dr]'\q\l + U2{^, t)) . 

Using these expressions as well as equations (4.18), equations (4.26) become 

Au = iUi, A12 = -iQn, A21 = A22 = -iU2. (4.30) 

Denoting by (^'i, ^'2)^ the second column vector of replacing A in equation (4.28) by equations 
(4.30), and noting that 72 = —7, where 7 is defined by equation (3.3), equation (4.28) implies that 
*i,*2 satisfy precisely equations (3.4) where the vector v is defined by equation (3.5). 
The vector (^1, ^'2)^ also satisfies 

= -^g*2, ^2^ = \q-^i- (4.31) 

It can be verified that the compatibility of equations (3.4) and (4.31) yields equations (1.3). 

It remains to: (a) show that the Fourier transform of equation (4.11) yields equation (1.8a); (b) 
justify the choice of the constants of integration of equations (4.22); (c) establish that q{0,rf,t) = 
q{^, 0, t) = 0. These interelated facts can be proved by "reversing" the relevant arguments used in 
section 3. 

QED 
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5 The t-part of the Lax pair for the non- homogeneous case. 

For simplicity of notation throughout this section we suppress the t-dependence. 

Proposition 5.1 Let the vector ^ satisfy equations (3.1) where q satisfies equations (1.3) on the 
quarter plane. Then the function tpiiVi k) = ^'1(0,77, k) satisfies the following equation: 



^'^ ^ 2 8 Jo ^ 



+ UlTpl 



/oo 
dl [7(fe 5{k -l,k)+ e{k - I, k)] ^i{r], I) 

-00 

-| d?7'#i + dv'q^i^i H\q\^ £' dfjq^i = 0, (5.1) 

where: q, q^ in the first bracket arc evaluated at (0, rj); \q\, qipi, q^ipi in the integrals with respect to 
drj' are evaluated at {0,r]'); qipi in the integral with respect to dfj is evaluated at (0,?}); 7 is defined 
by equation (3.3); and S,e are defined by 



00 roo 



Sil, k) = — / d^d^ g^(e, 0)q{^ + ^, 0) - q{^, 0)qr,{^ + ^, 0) 

oTT Jo Jo ^ 



e{l,k) 



327r 



roo roo roo 

/ / / d^did^'q{^,OMi + ^,0\'q{^ + i + e,Oy 
Jo Jo Jo 



-iki-ik^'-il^ 



(5.2) 



(5.3) 



Proof We first show that ^' satisfies equations (3.4) where v is defined by equation (3.5) with the 
additional term 

/oo 
dl {6{k -l,k)+ e{k - I, k)) *(^, rj, I). 
-00 

Indeed, proceeding as in the proof of Proposition 3.1 we find that the forcing of the equation 
satisfied by V2 involves the following additional terms 



roo roo 

qv{C,0) d^'qi^',0)e-"'^ -q{C,0) <g,(^', 0)e-^'=« 

Q /"OO fOO 

-3^^-(?,0) J di'\q{i'M' j, d^'W,0)e-^'^". 



(5.4) 
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For the computation of these additional terms we have used the following expressions, where for 
convenience of notation we have suppressed the fc-dependence: 

V'i^(e,0) = -ig(?,0)e-'=«, 

(5.5) 

^2r,r, 0) = kU^, 0)^1 , 0) + 0)^-1, (4, 0). 

Denoting the first bracket in (5.4) by 

POO roo 

/ deA{^,e)e~''^ = / dl6{k-l,k)e-''^, 

J( J-oo 

replacing k — I hy I and — ^ by ^, we find 

Jo J-oo 

which yields (5.2). Similarly denoting the second term in (5.4) by 

poo roo POO 

/ <A(e,0 / di"B{i")e-'^^ = dleik-l,k)e-'^^, 

Ji J^' J-oo 

replacing A; — I by I, ^" — ^' by ^, and ^' — ^ by ^, we find 

POO fOO ^ ^ ^ ^ POO 

/ d^A{^, ^ + d^B{^ + C + Oe-^^«+«) = / dle{l, k)e''^, 

Jo Jo J-oo 

which yields (5.3). 
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Equation (5.1) follows by evaluating equation (3.4a) at ^ = 0. For this computation we use the 
following equations, where again for convenience of notation we have suppressed the ^-dependence: 



^-2(0,77) = 1 + ySdri'q{0,v')MO,v'), ^^2,(0,7?) = ^g(0, 7?)Vi(0, r?), 



^i,,(0,r?) = -i9^(0, 77)^2(0, 77) - lq{0,v)4'2,{0,ri), 



V^i,^(0,r?) = -ig^ (0,77)^2(0,7?) - i|g(0,7?)|2V;i(0,77). 



QED 



Proposition 5.2 Let S{k,l,t) be defined by equation (2.4), where = ^'2 exp [— i/c^] , ^ is 
defined in terms of q in equations (3.1) and q evolves in time according to equation (1.3). Let S be 
defined in terms of S by equation (3.9). Then S satisfies the equation 

iSt + % + Sr,r, + {Ul + U2)S + d|Fi(^, |, t)S{l 7/, t) + /" ^^^2(7/, fj, t)Si^, 7?, t) = 0, (5.7) 

Jo Jo 



where 



m,i,t) = \ [%{i,0)qi^,0) - g(l,0)(?,(C,0)] - ^g(^,0)g(|,0)^/dC'|g(^',0)|2, 

^2(77,77,*) = ^fe(0,^)9(0, 77) -^(0,^)9^(0,7?)] -^9(0,77)9(0, 77)^'' d77'|g(0, 7701'. (5-8) 
Furthermore, S satisfies the boundary conditions 



S{^, 0, t) = Trq{^, 0), S{0, ri, t) = Trq{0, 77). 



(5.9) 



Proof The analogue of equation (3.12) is now 



f + iISdv'\q{0,v')\' „/l 



A;2 



(5.10) 
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Equation (3.11) yields 



/oo "I roo roo 

-OO J —OO Jo 



Furthermore, 



/OO PCX 
dk^i{0,r],t,k) = / 
-OO J—c 



dk 



since ^'i(0, 77, t, /c) is analytic for Imk < 0. 

Taking the Fourier transform of equation (5.1) and using the estimate (5.10) we find equation 
(5.7). In this respect we note that if we denote the bracket in the rhs of equation (5.2) by A(^,^), 
then the contribution of the term involving 7 is given by 



/OO /"OO 1 /"OO /"OO /"OO /"OO 

di dk-f{k-i,t)Miy^^ = ^ di dk d^' diAi^',i)t 
-00 J-00 oTT J-00 J -00 Jo Jo 



-iki-i(k-l)i'+ik£, 



Ml) 



where = "01 0- Similarly the contribution of the term involving S is given by 



"I roo POO POO POO POO 

- / dkdl / / d|«'9(l,0)|g(C + |,0|2g(| + ^ + ^^0)e-^'=«-^'=^-^(*^-'^^^^^^ 

61 J -OO J -OO Jo Jo Jo 



g(?,o) 
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POO POO 

/ di / di9{i - I - i)q{l 0)|g(| + 1, 0)|25(|, r,) 
Jo Jo 
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/ di [ diq{i,OMi + i,0)\^S{ir^). 
Jo Jo 



The area of integration is depicted in Figure 5.1a. Making the change of variables ^' = ^ + ^, ^ = 





Figure 5.1 
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the area of integration is mapped to the area depicted in Figure 5.1b. Thus the relevant integral 
becomes 



Also regarding the last term of equation (5.1) we note that taking its Fourier transform we find 



Changing the order of the integration we find that the relevant contribution involves the last term 
of F2, see equation (5.8). 



6 Discussion 

In the last decade considerable progress has been made in the understanding of boundary-value 
problems for integrable nonlinear evolution PDEs in one spatial dimension. For example, two of the 
articles in this special issue are concerned with this development [9], [10]. Here, to our knowledge 
for the first time, a boundary- value problem is solved in two spatial dimensions. 

It has been emphasised by the author that an important difference between initial- and boundary- 
value problems is the following: For initial-value problems one needs to perform the spectral anal- 
ysis of the t-independent part of the Lax pair only, while for boundary- value problems the spectral 
analysis of both parts of the Lax pair is needed. In this respect we recall the following developments: 

(i) It is interesting that the spectral analysis of the t-part of the Lax pair (in addition to that of 
the t-independent part), was first considered not for an equation in one spatial dimension but 
for an equation in two spatial dimensions, namely the DSI equation: In [4] equations (1.3) 
were solved on the plane, but with nontrivial boundary conditions at infinity. 



-g(C,o) 

16 Jo 





QED 



The authors of [4], rather than performing the explicit spectral analysis of the relevant t- 
dependent eigenvalue equation they made use of a certain completeness relation (of course 
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the derivation of this relation is based on the spectral analysis of the Schrodinger eigenvalue 
equation). 

(ii) In [11] the independent spectral analysis of the two parts of the Lax pair of the nonlinear 
Schrodinger equation led to the formulation of the solution in terms of two Riemann-Hilbert 
(RH) problems, which had to be solved in sequence. 

(iii) In [12] the above two RH problems were combined and the solution q was expressed in terms 
of a single RH problem. This RH problem has the distinctive and very useful feature of 
involving jump matrices with explicit exponential (x, t) dependence. 

(iv) It was shown in [13] that the above RH problem can be derived in a straightforward manner 
by performing the simultaneous spectral analysis of the Lax pair. The rigorous proof that 
the solution of this RH problem yields the unique solution of the given initial-boundary value 
problem was presented in [7], [8]. 

In the present article we have implemented for equations (1.3)-(1.5) the construction of (i) 
above. The main reason for using (i) instead of (iv) is the fact that for the case of homogeneous 
Dirichlet boundary conditions the analysis of the t-part of the Lax pair is very similar to the analysis 
presented in [4]. 

The implementation of (iv) for the initial-boundary value problem formulated in Proposition 
1.1 remains open. Furthermore, the analysis of the associated global relation, which characterises 
gi and /i in terms of go, fo and qo, also remains open. 
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